Spectral diagrams of Hofstadter type for the Bloch electron in three dimensions 
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Flux-energy and angle-energy diagrams for an exact three-dimensional Hamiltonian of the Bloch 
electron in a uniform magnetic field are analyzed. The dependence of the structure of the diagrams 
on the direction of the field, the geometry of the Bravais lattice and the number of atoms in an 
elementary cell is considered. Numerical evidence is given that the angle-energy diagram may have 
a fractal structure even in the case of a cubic lattice. It is shown that neglecting coupling of Landau 
bands changes considerably the shape of the diagrams. 
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The spectral properties of an electron in a two- 
dimensional (2D) periodic structure with the Bra- 
vais lattice A in the presence of a uniform magnetic 
field B are determined by commensurability or non- 
commcnsurability of two geometric parameters: the area 
S of an elementary cell of A and the square of the mag- 
netic length l\j = h/ /no (here \i is the mass of the elec- 
tron and u> is the cyclotron frequency). If S/2tt1 2 4 is 
a rational number, then the electron energy spectrum 
has band structure, otherwise the spectrum is a fractal 
set. As a result, the flux-energy diagram for the spec- 
trum has a remarkable recursive structure predicted by 
M. Ya. Azbel' and numerically discovered by D. R. Hofs- 
tadter in the framework of the tight-binding approxima- 
tion 0. The essential point in the appearance of such 
a structure is the size quantization of the electron mo- 
tion along the field direction; hence it seems likely that 
in three-dimensional (3D) periodic systems, the fractal 
structure of the flux-energy diagram must disappear 
Surprisingly, it was shown recently in a series of papers 
that a fractal structure is visible in the diagram depict- 
ing the dependence of the spectrum on the angle between 
the field B (with fixed strength B) and a fixed vector in 
A, if A is an anisotropic rectangular lattice. Moreover, 
a series of energy gaps as in Hofstadter's butterfly arise 
in the isotropic case unless B points in high-symmetry 
crystallographic directions 0. 

The tight-binding approximation method used in 0, 
0, is based on a series of considerable simplifica- 
tions of the initial periodic Landau operator, which im- 
poses severe limitations of the method. In particular, 
the models considered in and give no way to take 
into account the effects of the interaction between Lan- 
dau levels, which have profound effects on the shape of 
flux-energy diagrams and on the accompanying integer 
quantum Hall conductivity 0. Even without the con- 
sideration of the Landau band coupling, the flux-energy 
diagrams obtained for the periodic Landau operator dif- 
fer considerably from those for the tight-binding model 
6]. In this paper, we get rid of the restrictions imposed 
by the tight-binding approximation method and analyze 



the angle-energy and the flux-energy diagrams for the 
3D Landau operator perturbed by periodic potentials of 
various geometries. To deal with an explicitly solvable 
model, we consider periodic perturbations represented as 
sums of short-range potentials; in the zero-range limit 
we get an explicit expression for the dispersion relations, 
which is very useful for the numerical analysis. Such po- 
tentials were already used in the context of spectral prob- 
lems for a 3D magneto-Bloch electron , ■ In the 2D 
case, the flux-energy diagrams for zero-range potentials 
were considered in |9j- 

The Hamiltonian H being studied is the sum H = 
H + V, where 
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is the Landau operator with the vector potential A(r) = 
B x r/2, and V is the potential of a crystal lattice T with 
the Bravais lattice A. Let K be the set of all nodes of T 
lying in an elementary cell of A, then V is represented in 
the form 



nr)^^K(r-A). 
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where V K (r) is the confinement potential of the node k 
which is supposed to be short-range. More precisely, we 
choose V K (r) = c K W(r — k), where W(r) ~ outside a 
small sphere of radius R centered at zero, J W(r) dr = 1, 
and the coupling constant c K is of order R. At the zero- 
range limit, R — > 0, the potential V K is characterized 
by one parameter only, namely, by the scattering length 
p K , which is related to the binding energy of the ground 
state for V K by E K = —h 2 /2fj,p 2 l 10]. Moreover, at this 
limit the Green function G of H, G(r,r';E) = (r\(E - 
-ff) -1 |r'), has the following explicit expression in terms 
of the Green function G (r,r';E) of H [l^, 0: 

G(v,t , ;E) = G (t,v , ;E)- £ G (r, T ,E) 

7, -e'er 

x(Q-\E)) G ( 7 ',r';E). (3) 
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Here Q l (E) is the matrix inverse to the infinite matrix 
Q(E) with elements 
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+(l-5 7i7 /)G (7,7';-E), (4) 



where GJj 011 denotes the renormalization of the unper- 
turbed Green function Go- 



Gr(r,r';E)^G (r,r';E) 



p cxp [— i7rb(r x r')] 
2ith 2 |r-r'| 



with b = Be/27rftc (note that <po = 2-Khc/e is the mag- 
netic flux quantum, therefore b is the density of the mag- 
netic flux). The explicit form of Go is well known [llj . 
El: 



G {r,r';E) = $(r,r')F(r-r';E), 
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where 
<P(r,r') = 

F(r;E) = 
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Here ru is the projection of r on the direction of the 
field B, = r — rn, Li(x) denotes the l-th Laguerre 
polynomial, ande ; = Hu(l+l/2). Note that G cn (r, r'; E) 
is well defined at r = r' and independent of r: 
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where C( s i2ii_i s t ne generalized Riemann (or Hurwitz) 
£- function |l3| . 

It is clear that at least for E < fku/2 the energy E 
belongs to the spectrum of H if the matrix Q(E) is not 
invertible. Using the irreducible representations of the 
magnetic translation group (MTG) for H the prob- 
lem to invert the infinite matrix Q(E) may be reduced to 
a problem of finite-dimensional algebra. This reduction 
requires the so-called " rationality condition" : the field 
B is said to be rational with respect to the lattice A, if 
for a basis ai, a2, a3 of A, the numbers h(&j x a&) are 
rational ^ij. In this case a basis a.j (j = 1,2,3) can 
be chosen in such a way that 77 = b(ai x a2) > and 
b(a2 x a 3 ) = b(ai x a 3 ) = 0. Let n — N/M, where 
N and M are coprime positive integers. Then all the 



irreducible representations of the MTG which are triv- 
ial on the center of the group are M-dimensional and 
are parameterized by a 3D torus T 14]. It is conve- 
nient to choose coordinates of a point p from T such 
that < pi < M -1 , < P2,P3 < 1- Then we can form 
the following (MK) x (MK) matrix Q, where K is the 
number of nodes in K: 

Q q , q '(p, E) = exp[-7riTO'b(K' X Sl 2 )] 

00 

x V" exp \-m{2\ ■ p + r]X 1 (MX 2 + to))] 

Ai .A2 , A3 — — oc 

xQ(Aiai + (X 2 M + TO,)a 2 + k, m'a 2 + k' , X 3 a 3 ;E) 

x exp [7ri(Aiai + (X 2 M + m)a 2 )(b x k)] . (8) 

Here q denotes the pair (m, k) with k G K, to = 
0, . . . , M — 1. Now the dispersion relation for H reads 



detQ(p,£) = 0. 



(9) 



Equation @ has for fixed p infinitely many solutions 
E s (p) (dispersion laws), which are continuous with re- 
spect to p; each eigenvalue E s (p) is M-fold degenerate. 
By definition, the magnetic miniband J s is the set of all 
values of E s (p). The minibands J s lying below Huj/2 
form a piece of the spectrum of H which may be at- 
tributed to broadening the ground state of Ho; therefore, 
this piece is nothing but the lowest Landau band. Ac- 
cording to Eqs. © and ©, this band consists of KM 
magnetic minibands which can overlap. If the overlap- 
ping is absent, then the lowest Landau band approaches 
a Cantor set as rj approaches an irrational number (and, 
therefore, M tends to infinity). 

Using Eqs. JSJ and © we analyze numerically the 
structure of the lowest Landau band for various types of 
crystal lattices and ranges of the field B employing two 
ways to force n to approach an irrational number: (1) we 
change the value of B keeping the direction of B fixed; 
(2) we change the direction of B keeping the value of B 
fixed. The lattice constant a — 7.5 nm is chosen relevant 
to the geometric parameters for the 3D regimented quan- 
tum dot superlattice considered recently in ^5|. As to 
the scattering length, we put p ~ 1 nm, this corresponds 
to the binding energy E ~ 30 meV. 

The flux-energy diagrams for the simple-cubic lattice 
under various directions of B are depicted in Fig. 1 (here 
p = 2.5 nm). If the magnetic field is directed along an 
edge of the cubic elementary cell, then all subbands over- 
lap and the gaps are absent (Fig. la). If the symmetry of 
the magnetic field with respect to the lattice decreases, 
then the number of open gaps increases and has ten- 
dency to infinite magnification (Figs. lb,c) in full agree- 
ment with 0]. Fig. 2 shows the angle-energy diagrams 
for the lattice with the same geometric parameters and 
p = 1.1 nm (in Figs. 2-4 77 = 1; for a = 7.5 nm this 
corresponds to the field strength ~ 70 T). If the field B 
is rotated inside a face of a cubic elementary cell, then 



1/4 1/2 3/4 1 
1 




(b) 



-0.5 



■1.0 



W -1.5 



-2.0 



-2.5 



! 



—i i i i_ 



1/5 1/3 1/2 2/3 4/5 
11 



(C) 



-0.5 



-1.0 



69 

W -1.5 



-2.0 



-2.5 



I 'l'-lfn , , 



i i iiiiii ii in ] 



| jilil 



lit" 



f IE 







_i i_ 



1/5 1/3 1/2 2/3 4/5 

11 



15 30 45 60 75 90 



(b) 



-0.5 



W -1.0 



-1.5 



|!jll|!!|;iii 



illllll I 



:; l ! ill 



(C) 



2.0 



-4.0 



-6.0 



—I I I l_ 



15 30 45 60 75 

e 



_i i i i i i i 

15 30 45 60 75 90 



90 



FIG. 1: The lowest Landau band of the simple-cubic lat- 
tice plotted against the magnetic flux for various directions 
of B: (a) B = (0,0, B); (b) B = (|B,0,|B); (c) B = 
(i§B, £gB, f B). Here and below E = h 2 /^a 2 . 



FIG. 2: The lowest Landau band of the simple-cubic lat- 
tice plotted against the tilting angle 8 for various tilting 
orientations of B: (a) B = 73(sin 8, 0, cos 8); (b,c) B = 
B(| sin 6, | sin 8, cos &), in panel (c) coupling between Lan- 
dau levels is neglected. 



the gaps, in general, overlap and the diagram reveals no 
fractal structure (Fig. 2a). On the other hand, if the ro- 
tation plane forms a dihedral angle (p — atan(3/4) with 
the plane of the face, then the fractal structure of the 
diagram is clearly visible (Fig. 2b). Therefore, the con- 
dition of anisotropy 0] is not necessary for the appear- 
ance of a fractal structure in the angle-energy diagram. 
This can be understood taking the limit B — > oo in the 
matrix Q. Eq. (jSJ shows that in this limit Q approxi- 
mates the matrix of the tight-binding Hamiltonian from 
3] with energy- depending coefficients tj, and this depen- 
dence leads to an anisotropy of the limiting matrix with- 
out any anisotropy of the crystal lattice. 

If in the representation |JBJ we restrict ourselves only to 
the first term, then this projection of the Green function 
on the lowest Landau level models the system without 
interaction between Landau levels. The corresponding 
diagram is given in Fig. 2c. 



To compare our results to those from |3j we con- 
sider the angle-energy diagrams for the tetragonal 
monoatomic lattice with p = 2.5 nm (Fig. 3). If the 
field B is rotated in a face of an rectangular elementary 
cell, then we see a typical ID-like energy spectrum in full 
agreement with It is interesting that the chemical 
anisotropy radically transforms the shape of the angle- 
energy diagram. To show this, we consider a tetragonal 
double-atomic lattice with two sorts of atoms in an el- 
ementary cell with scattering lengths p\ — 2.5 nm and 
P2 = 0.6 nm (Fig. 4). If the magnetic field is rotated in a 
face of anisotropy, then the angle-energy diagram looks 
like deformed Hofstadter's butterfly. Since in this case 
K = 2, there is a doubling of minibands that causes the 
appearance of a wide gap in the diagram (" atomic" gap) . 
On the other hand, if B rotates in the plane perpendicu- 
lar to the anisotropy axis, only the atomic gap appears in 
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field (Eqs. iJSJ and (JSJ)) and analyzed the corresponding 
flux-energy and angle-energy diagrams. For anisotropic 
rectangular crystals we confirm the results of 0] and 
concerning the structure of these diagrams. Moreover, 
we show that the anisotropy is not necessary for the ap- 
pearance of a multitude of gaps both in flux-energy and 
angle-energy diagrams. Neglecting coupling of Landau 
bands leads to a substantial deformation of angle-energy 
diagram. Chemical anisotropy of the considered crystal 
(the presence of distinct sorts of atoms in an elementary 
cell of the Bravais lattice) leads to another deformation 
of the diagram, in particular, to the appearance of an 
additional wide gap. 

The work is partially supported by DFG SFB-288 "Dif- 
ferentialgeometrie und Quantenphysik" and by Grants of 
DFG-RAS, INTAS, and RFBR. 



FIG. 3: The lowest Landau band of the tetragonal 
monoatomic lattice (with lattice constants (a/2, a, a)) plot- 
ted against the tilting angle 9 for various tilting orientations 
of B: (a) B = B(sin 9, 0, cos 9); (b) B = B(0, sin 9, cos 9). 
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FIG. 4: The lowest Landau band of the tetragonal double- 
atomic lattice (with lattice constants (2a, a, a)) plotted 
against the tilting angle 9 for various tilting orientations of 
B: (a) B = B(0, sin 9, cos 9); (b) B = B(sin 9, 0, cos 9). 



the diagram: (Fig. 4a arises from Fig. 2a by the doubling 
of all bands). 

As a conclusion, we have derived a new dispersion re- 
lation for a 3D Bloch electron in a uniform magnetic 
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